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Abstract. We propose a procedure to evaluate the impact of nonlinear couplings on the 
evolution of massive neutrino streams in the context of large-scale structure growth. Such 
streams can be described by general nonlinear conservation equations, derived from a multiple- 
flow perspective, which generalize the conservation equations of non-relativistic pressureless 
fluids. The relevance of the nonlinear couplings is quantified with the help of the eikonal 
approximation applied to the subhorizon limit of this system. It highlights the role played by 
the relative displacements of different cosmic streams and it specifies, for each flow, the spatial 
scales at which the growth of structure is affected by nonlinear couplings. We found that, at 
redshift zero, such couplings can be significant for wavenumbers as small as k = 0.2/i/Mpc 
for most of the neutrino streams. 
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1 Introduction 

Statistical properties of the large-scale structure of the universe have long been proposed 
as an efficient instrument to constrain cosmological parameters. In this context, a careful 
account of the role played by massive neutrinos is crucial. So far, it has often been overlooked 
in the nonlinear or quasilinear regime because of the technical complexity specific to the 
study of massive neutrinos (see recent attempts in [1-3]). This is all the more unfortunate 
that cosmological observations can fruitfully improve our knowledge of those particles. The 
signature of neutrino masses on cosmological observables is indeed expected to be significant 
enough for those masses to be constrained observationally [4-9]. 

In the linear regime, the effect of neutrinos is now well understood (see refs. [10-12]). 
The need for nonlinear corrections in their equations of motion has been raised because the 
cosmological observations that are the most sensitive to neutrinos masses, i.e. for wavenum¬ 
bers in the 0.1-0.2 h/Mpc range, precisely correspond to the mildly nonlinear regime. To 
deal with this issue, several strategies can be adopted. Until recently, in analytic works, the 
neutrino fluids had always been treated in the linear regime, nonlinear couplings being intro¬ 
duced in the dark matter description only. This nonlinear treatment can be implemented with 
the help of the Renormalization Group time-flow approach [1, 13]. Improvements upon such 
schemes have been proposed in [3], which consists in a hybrid approach that matches the full 
Boltzmann hierarchy to an effective two-fluid description at an intermediate redshift. [14] is, 
for its part, a systematic perturbative expansion of the Vlasov equation in which high-order 
corrections to the neutrino density contrast are computed without the explicit need to track 
the perturbed neutrino momentum distribution. 

Ideally, however, the fully nonlinear evolution of the neutrino fluid should be depicted. 
A natural way to do so would be to take inspiration of the standard linear description, which 
relies on the Boltzmann equation, and extend the harmonic decomposition of the phase-space 
distribution function to the nonlinear regime. This has been done in [15] but this method 
turned out to be particularly difficult to handle. In [16, 17], we proposed to describe massive 
neutrinos as a superposition of single-flow fluids, the equations of motion of each of them 
being written in the nonlinear regime. 

In this paper, we are interested in exploiting those theoretical developments in order 
to identify the scales at which nonlinear couplings in the neutrino fluids are expected to 
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play a significant role. In order to do so, we apply the eikonal approximation to the nonlinear 
eqnations of motion. Note that this approximation had already been exploited in the literature 
to develop a Perturbation Theory approach for cold dark matter ([18]) and had proved to be 
able to capture the leading coupling effects. 

The article is structured as follows. In section 2, we recall the form of the nonlinear 
equations of motion describing the time evolution of non-interacting fluids, relativistic or not, 
when using a multi-fluid approach. Section 3 explains in detail how the eikonal approximation 
can be implemented in those equations and emphasizes the key role of relative displacement 
fields. Finally, in section 4, power spectra of the relative displacements between neutrino 
fluids and the cold dark matter component are presented. The impact on the growth of 
large-scale structure is then discussed in a quantitative way. 


2 Nonlinear equations of motion (multi-fluid description) 


Following [16, 17], it is now clear that any non-interacting relativistic fluid can be divided 
into several flows, each of them evolving then independently until first shell-crossings. In 
cosmology, this approach obviously applies to massive neutrinos since they are free-streaming. 
In this framework, each flow can be defined as the collection of all the particles (with a mass 
m) having the same comoving momentum. They are entirely characterized by two coupled 
fields, namely the comoving number density Uc and momentum P^. Those fields obey the 
following general equations. 


d d {\ 


= 0 , 


( 2 . 1 ) 


with P>^{r],x^) = Pi,{rj,x^) and P^P^ = —m?, being the metric, and 




( 2 . 2 ) 


These equations directly ensue from the matter and momentum conservation equations, ap¬ 
plied to each flow. At this stage, no perturbative expansion of the metric is involved. The 
properties of the whole fluid are then inferred by examining an appropriate number of such 
flows, each of them being labeled by the initial value of its field Pi, denoted (found to be 
constant at zeroth order in Perturbation Theory). The initial number density of particles in 
each flow is constrained by the choice of initial conditions. For instance, the case of adiabatic 
initial conditions is described in detail in [16]. 

We are interested here in equations i) involving linearized metric perturbations (but 
non-linearized fields) and ii) rid of the coupling terms that are subdominant at subhorizon 
scales. They can be written in terms of the number density contrast 
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and of the velocity divergence field (in units of —H, % being the conformal Hubble constant) 
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since the field Pi has found to be potential in this regime^. For a generic perturbed Friedmann- 
Lemaitre metric^ , whose time variable is the conformal time rj, 

ds^ = a^{r]) [—(1 + 2A)dr]^ + 2Bidx’'dr] + {6ij + hij)dx^dx^~\ , (2-5) 


one has in Fourier space for the mode characterized by the wave vector k ([17]) 
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Besides, 5T-j(k) is a source term given by 
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These equations contain also the generalized kernel functions, adapted to relativistic flows. 


ai?(ki, k2; r) = (5Dirac(k - ki - k2) 
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As mentioned in [17], they are extensions of the kernel functions found for pressureless fluids 
of non-relativistic species (see [19] for details in this context.). 

In any practical implementation, it is necessary to consider a collection of N streams. 
In that case, the general equation of motion can conveniently be written in terms of the 
time-dependent 2A^-uplet, 


T,(k) = (<5.,(k),0.,(k),..., ,5.^(k),0.^(k))^. (2.12) 


Before shell-crossing, it can incorporate all the relevant species (neutrinos, dark matter, 
baryons) as long as they interact only via gravitation. In this context, the equations (2.6) 
and (2.7) of all the flows can formally be recast in the form^ 

5,T,(k) + Tb(k) = 7„''^(ki, k2)T6(ki)Te(k2), (2.13) 

^This property, rigorously demonstrated in [17], generalizes that of non-relativistic species. 

^Units are chosen so that the speed of light in vacuum is equal to unity. 

^The Einstein notation for the summation over repeated indices is adopted and, in the right hand side of 
this equation, it is assumed that the wave modes are integrated over. 
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where djdz = adjda and the indices a and b run from 1 to 2N. The matrix elements 
gather the linear couplings. They contain in particular the way in which the source terms 
5T-j(k) can be re-expressed as a function of the 2A^-uplet elements. The left hand side of this 
equation is nothing but the linear held evolution. The right hand side contains the coupling 
terms. More precisely, the symmetrized vertex matrix 7 „^'^(ki,k 2 ) describes the nonlinear 
interactions between the Fourier modes. It is given by 


with 7 „^''(ki,k 2 ) 


72p-i ^(ki,k2)- a/j(ki,k2, Tp), 

2(Tp)o 

(2.14) 

72/^(ki,k2)- , /3R(ki,k2, Tp), 

v'T’lo 

(2.15) 

(k 2 ,ki) and 7 ^^'^ = 0 otherwise. 



3 The eikonal approximation 

The eikonal approximation^, developed in [18], is based on the observation that the amplitudes 
of the kernel functions describing mode couplings, an and Pr, signihcantly depend on the 
ratio between the wave numbers at play. This has been observed for non-relativistic huids 
and we show here that it is the case for relativistic huids too. This property leads to the 
idea that the right hand side of eq. (2.13) can be split into two integration domains. One is 
called the hard domain and encompasses modes whose wavelengths are of the same order (and 
for which the coupling functions are always hnite). The other one is referred to as the soft 
domain. It is made of modes of very different wavelengths for which the coupling functions, 
of the order of the wavelength ratio, are large. 

The main idea is that there are regimes in which the dominant coupling structure is 
in the soft domain. This is the case for instance in early-time fluids containing baryons and 
cold dark matter. At the time of recombination, the baryon velocity drops steeply whereas 
the velocity of cold dark matter is not affected by decoupling. It means that, at intermediate 
scales (i.e. between the Silk damping length and the sound horizon), the relative velocity Uj-ei 
between the baryon flow and the one of cold dark matter is substantial. Because of this, for 
scales at which k > aH /(at decoupling), the wavelength of gravitational potential 
wells is too small for baryons to fall in them before being pushed towards another direction. 
Eventually, this phenomenon induces a damping of the matter power spectrum. The relative 
motion between cold dark matter and baryons and its effects on the matter power spectrum 
have been highlighted in [21, 22]. Such studies illustrate the relevance of couplings between 
large scale modes and small scale modes in the framework of structure formation. For instance, 
in the case of mixtures of baryons and cold dark matter, the typical coherence length of the 
relative velocity held is of the order of few Mpc, which is much larger than the scales at which 
basic baryonic objects start to form under gravitational clustering (of the order of 10 kpc). 
The same formalism can be used to obtain the large k behavior of the propagators in case of 
a single pressureless huid, reproducing the results obtained in [23, 24]. In the same spirit, we 
propose here to investigate the impact of the relative motion between given neutrino streams 
and the cold dark matter huid with the help the eikonal approximation. It will allow us to 
infer the amplitude of neutrino coupling effects. 

^In this context, the term refers to diagram resummations performed in quantum electrodynamic field 
equations, [20]. 
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For convenience, let us assume that the soft domain is obtained for k 2 <C ki in eq. (2.13). 
We have then k = ki and the contribution corresponding to the soft domain can be viewed 


as a mere corrective term in the linear equation describing the evolution of the mode k. In 
other words, eq. (2.13) can be rewritten 


k) + nj’iz, k)^'fe(2;, k) - -^(z, k)^'b(z, k) 

= 7a'"(ki,k2) ^ 6 (ki,^) ^,{k2,z)] (3.1) 

-I rC 


with 


z)^2j^ d^q “"■ 7 .''=(k, q)'k,(q, z) . (3.2) 

The soft momenta q (i.e. g <C A:) at play in eq. (3.2) are integrated over so that Ej^{z,k) is 
independent on Ta(z,/c). It is a mere time and scale dependent matrix. The fact that the 


integration domain is restricted to the soft wave numbers in eq. (3.2) is the key element. 


Conversely, in the right-hand side of eq. (3.1), the implicit convolution product excludes 


the soft domain (i.e. all the modes concerned have comparable wavelengths). When the 


contribution of the hard domain is negligible, eq. (3.1) ean then be viewed as the equation 


of motion of the mode k evolving in a medium perturbed by large-seale modes. It therefore 
encodes the way in which long-wave modes alter the growth of structure. Once rf^{z,k) is 
given, eq. (3.1) can be solved as a linear equation. This is precisely the eikonal approximation 
of the global equation of motion. 


In practice, applying the eikonal approximation to Ej’[k, z) means that the vertex values 
that appear in this quantity have to be computed assuming k 2 ki. In this framework, one 
deduces from eqs. (2.10), (2.11), (2.14) and (2.15) that the eikonal limit of the vertex elements 
is 



-S 2 'p(k,k 2 ) = -h 2 V 2 , TT^^k- (k 2 -^^r;) , (3.3) 



(3.4) 


This expression depends on each flow through its initial momentum 77 , which vanishes in 
the standard non-relativistic equations. Besides, (Tp)o —>■ —ma in the non-relativistic limit so 
that one recovers the expected formulae in this limit (see [21, 25]). 

In the following, we exploit the consequences of this approximation in order to evaluate 
the relevance of the coupling terms. As a hrst step, we can notice from eqs. (3.2), (3.3) and 

(3.4) that the two non-zero elements coming from the flow labeled by Tp in the H matrix, 
E 2 p^ and E 2 pf(f^, are proportional to the velocity divergence of large-scale modes. Thus we 


can write 



\k.dp{zQ,z) h, 


b 


(3.5) 


where a and b are either 2p or 2 p — 1 and where dp is the total displacement held induced by 
the large-scale modes in the huid labeled by Tp. It reads necessarily 
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Note that this displacement is superimposed on the zeroth order displacement field induced 
by the homogeneous momentum of the flow. This background displacement field is given by 
(see eqs. (2.6) and (2.7)) 



We will discuss this point in greater detail in the next section. 

The impact of the eikonal correction introduced in the equation of motion depends on 
the way in which the various large-scale modes contribute to the displacement fields dp. 
The nature of the displacement is crucial in this context. In particular, one expects global 
displacement fields, which affect all species in a similar way, to induce a mere phase shift in the 
solution of the eikonal limit of (3.1). Remarkably, such displacements have no impact on power 
spectra, provided that the fields at play are evaluated at the same time (see the end of this 
section). Using the standard language of cosmology, we call those particular displacements 
“adiabatic displacements” in the following. On the contrary, relative displacements between 
species can induce a damping in power spectra. It is simply due to the fact that species 
must evolve in phase (at least during a small period of time) for couplings between them to 
generate a significant growth of perturbations. This phenomenon has been highlighted for 
the first time in [21] an reconsidered in [22]. In this paper, we show that such considerations 
can be extended to the study of relativistic fluids. 

Knowing this, it is convenient to decompose displacement fields into 

dp{z, zo) = dadiab(^, zo) 5dp{z, zo), (3.8) 

where dadiab denotes naturally the adiabatic part. The other part, 6dp, represents relative 
motions between fluids. For non-relativistic species, it is known that the most growing mode 
is part of the adiabatic modes. Since neutrinos become non-relativistic at late time, one 
expects the most growing mode of each fluid of neutrinos to become also an adiabatic mode 
ultimately. 

When considering only the adiabatic part of the displacement field, the solution of the 
eikonal equation of motion is easy to find. Its form is related to the extended Galilean invari¬ 
ance of the equations of motion, well established for non-relativistic species and uncovered for 
the system (2.6)-(2.7) in [17]. Indeed, it has been shown that this system is invariant under 
the following transformations. 


X* = X* -k di{z), (3.9) 

z = z, (3.10) 

^a{z,^) = (3-11) 

where the last transformation can equivalently be written '^^(r), k) = exp(ik • d(?7))^a(77, k). 
It means that a homogeneous time-dependent displacement which disrupts the medium can 
be re-absorbed in a global phase shift of the Fourier transforms of the fields. Here, we are 
interested in making the large-scale adiabatic displacement fields play the role of the disturbers 
of the medium. 

More explicitly, it is known that the solution of the standard linear system can be fully 
described with the help of its Green function, g^^(z,zojk), defined in such a way that (see 
appendix A for details) 

Ta(z,k) = 5f„^(z,zo;k)^'fe(zo,k), (3.12) 
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with z and zq two arbitrary times. Besides, when the displacement held is purely adiabatic, 
the Green functions of the corrected linear system (i.e. the linear system in which an 
eikonal correction has been added) are related to those of the naked theory by 

L’’iz,Zo;'k) = 5r„^(2;,zo;k)exp(ik.dadiab(2:,^o))- (3.13) 

A direct consequence of the symmetry property is that, after a transformation of the form 
(3.11), ensemble averages of any product of helds V’a(' 2 ^a)ka) become proportional to 
exp[i^^ka • dadiab('2a)]) which is unity when the helds are computed at equal time. Indeed, 
statistical homogeneity imposes Yla ka = 0- This is the reason why equal time spectra or 
poly-spectra are not sensitive to the presence of adiabatic displacements. 

Yet, we are interested in all contributions to the displacement held, including those 
that induce large motions between species. As already mentioned, in early-time mixtures of 
baryons and cold dark matter, those relative displacements (i.e. those non-adiabatic contri¬ 
butions that develop in the nonlinear regime) are the leading contributions to the nonlinear 
evolution of the power spectrum. To sketch the impact of massive neutrinos on the nonlin¬ 
ear growth rate, we evaluate in the next section the amplitude of the relative displacements 
involving them and we compute the corresponding power spectra. 


4 Relative displacements and power spectra: quantitative results 


As stressed in the previous section, in the eikonal approximation, the displacement helds 
that can cause a damping of the growth of structure are those that differ from one huid to 
another. Hence, in this section, we compute on the one hand the power spectrum of the 
total displacement held of cold dark matter and, on the other hand, the power spectra of the 
relative displacement helds of relativistic hows (with respect to the a priori dominant cold 
dark matter component). Note that, in perturbation theory, such displacements are sums of 
terms of different orders (the zeroth order contribution dp^^ being dehned in eq. (3.7)), 

dp(2:,^o) = d®(2;,2;o) djp\z,zo) + ... (4.1) 

Note also that the damping due to non-homogeneous corrections can be signihcant only if the 
considered how is non-relativistic (since in that case dr^^( 2 ;, zq) is small) or if k is orthogonal 
to the zeroth order contribution (since in that case k.dr*^^(z, zq) is small compared with the 
other contributions). Hence, in the following, we compute the expected values of k.(dp — dcdm) 
as a function of the angle between the initial momentum r of the considered huid and k. 
Calculations are performed in the conformal Newtonian gauge, 

ds^ = a? (rj) [— (1 -|- 2'0) drj^ + (1 — 2</)) dx'^dx^dij] . (4.2) 


This choice will allow us to take advantage of the numerical work presented in [16]. First, let 
us dehne two transfer functions Dcdm as 



ma 



dz 0cdm , q) 


'20 


F>^"i(2;,2:o,q)V’init(q), 

Tlcdm('Z) ^0) H)V’init(q)) 


(4.3) 

(4.4) 
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where V’init is the initial value of the potential ip. The statistical properties of those quantities 
are entirely encoded in their initial power spectra Pipiq), defined so that 

(V'init(q)V'init(q')) = (2vr)^(5Dirac(q + q') P'lpiq)- (4.5) 


Using the expressions of the transfer functions, the contribution of each mode to the 
displacement field of each relativistic flow reads 


,{z,zo-,q) = -i (^^D^^\z,zo,q) - {z, zo,q)^ V'init- (4.6) 


Besides, for the cold dark matter component, one simply has 


dcdm(’2^? ^ 0 : Q) 





(4.7) 


Furthermore, one can notice that the displacement field of a relativistic flow along an 
arbitrary direction k depends on the angles between both q and r and k and r. After 
integration over the other angles, one finds for the variances of respectively k.dcdm and k.(dp — 
dcdm)) 

((k.dcdm)^) = 47rfc^ J dqP^{q) ^ |T>cdm(q)|^ (4.8) 

and 


((k.(dp - dcdm))') 


dqP^iq) J ^d/i 


= 27Tk'^ X (4.9) 

^(1 - Mfc)(l - fJ-'^)\D^p\q) - T>cdm(q) - .^(^^(q)^ + /rVil-Di^^(q)P , 


where Hk is the Cosine of the angle between the initial momentum of the flow Tp and k and 
where an integration is made over /r. Cosine of the angle between Tp and q. In eq. (4.9), 
one can notice that the dependence of the r.m.s. with respect to /i^ is such that it does not 
vanish either for an initial momentum Tp orthogonal to k (i.e. when = 0) or for an initial 
flow momentum along k (i.e. when /x^ = 1). 

On Figure 1, we present the per mode contribution to the right hand side of eq. (4.9) 
in the particular case of neutrino fluids for /Xfc = 0 on the left panel and /Xfc = 1 on the right 
panel. The results have been computed assuming a single species of neutrinos whose mass is 
nil, = 0.3 eV and using the cosmological parameters derived from the Five-Year Wilkinson 
Microwave Anisotropy Probe (WMAP 5) observations. Besides, the values of the initial 
power spectra have been obtained under the assumption that, in the cosmological model we 
adopt, metric fluctuations are initially adiabatic and characterized by the scalar spectral index 
Us ~ 0.96. We can see that the resulting neutrino power spectra are comparable in amplitude 
to the cold dark matter one (represented by a thick dashed line). 

Denoting = 1/A:((k.dcdm)^)^'^^ and ad^ = 1/A:((k.(dp - dcdm))^)^'^^ we find that 

is of the order of 2 to about 5/i“^Mpc for the cosmological modes we used for which 
fjdcdm ~ 6/x“^Mpc,. This result depends of course on the flow considered and is actually of 
the order of the cold dark matter value. What does it mean? Similarly to what happens 
in mixtures of baryons and cold dark matter, one expects the perturbation growth to be 
damped for wave numbers larger than or comparable to 1 /ad ^, for which k • d^^^ is expected 
to be hnite. Besides, such a non-adiabatic damping is potentially larger than the damping 
due to the homogeneous displacements dp^^ of each flow when /x^ is close to zero. A precise 
determination of the amplitude of these effects would require a full analysis of the nonlinear 
evolution of the system. We leave this for a future study. 
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Figure 1. Power spectrum of the relative displacement as a function of the mode q for different 
neutrino flows. The quantities are defined so that the right hand sides of eqs. (4.8) and (4.9) 

read f dqP^^cdm(q)- The values of t range from 2.25 ksTo (bottom line) to 18 ksTo (top line). 

On the left panel /ifc is set to 0, on the right panel fXk is set to 1 and the neutrino mass is set to 0.3 eV. 
The gray dashed line represents the power spectrum of the cold dark matter displacement. 


5 Conclusion 

Describing neutrinos as a collection of single-stream fluids is an efficient strategy to infer 
their impact on the growth of the cosmic structure. Using this approach, one indeed gets a 
complete set of equations of motion that incorporate all the nonlinear effects of relativistic 
or non-relativistic particles. In the subhorizon limit, the system takes the form of eq. (2.13), 
which can be easily handled by a formalism originally developed to depict non-relativistic 
species. 

In this paper, we evaluated the amplitudes of the nonlinear couplings and determined 
for each flow the scales at which they are expected to impact significantly on the structure 
growth. For that purpose, we implemented the eikonal approximation into the general equa¬ 
tion of motion. We concluded that the impact of large-scale modes on an arbitrary mode are 
entirely driven by large-scale displacement fields whose expressions are given in eq. (3.6). The 
comparison between the displacement field associated with each flow of neutrinos and the one 
associated with cold dark matter makes easy the comparison between the power spectra of 
the relative displacements between neutrinos and cold dark matter and the power spectrum of 
the displacement of cold dark matter alone. We found as a preliminary result that couplings 
involving massive neutrinos (with a 0.3 eV mass) are expected to induce a damping of the 
perturbation growth in neutrino flows for wave numbers larger than (or of the order of) about 
0.2 to 0.5 /i/Mpc. A detailed quantitative analysis of the consequences of this phenomenon 
is yet to be done but those findings confirm the significance of nonlinear couplings in the dy¬ 
namical evolution of neutrino fluids. This sets the stage for further numerical studies beyond 
the linear regime. 

Acknowledgements: This work is partially supported by the grant ANR-12-BS05-0002 
of the French Agence Nationale de la Recherche. 
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A Integral form of the overall equation of motion 

For a finite number of flows, the overall equation of motion (2.13) can formally be written 
in an integral form. It requires the use of the associated Green operator, gj’{z, zq] k), of the 
linear system. This operator satisfies 

Ta(2;,k) = gJ’{z,zo]k)'i'h{zo,k), (A.l) 

with z and zq two arbitrary times. It is besides solution of the differential equation 

zq; k) + k) gj’iz, zo; k) = 0 (A.2) 

with the condition 

9aizo,zo;k) = 6j', (A.3) 

Sj’ being the identity matrix. Formally, the Green function is the ensemble of all the inde¬ 
pendent linear solutions of the system^. Denoting u'a\z,k) these solutions, g\ reads 

gl{z,zo,k) = ^4“)(2;,k)cf„)(zo,k), (A.4) 

a. 

where the variables c^^^(zo)k) are set so that (A.3) is satisfied. 

Studying in detail the Green operator of such a system is beyond the scope of this 
appendix. Suffice to note here that, unlike the case of a single pressureless flow, the Green 
operator generally depends on the wave mode k. This dependence is expected to gradually 
decay over time and to disappear at very late time, when all the flows have become non- 
relativistic. At this stage, the situation is then identical to the one of a collection of cold dark 
matter fluids. 

As for the standard system of non-relativistic particles, the knowledge of the Green 
operator of the equation of motion allows to write a formal solution (see [26-28]), which is 
given by 

T„(k,z) = gJ'{k,z,zo) T6(k,2o)-k 

+ [ dz' gj'{k,z,z') %‘'‘^(ki,k 2 )^ciki,z')^dik 2 ,z'), (A.5) 

JZQ 

with 'I'a(k, Zq) the initial conditions. Many of the approaches developed in order to improve 
upon standard Perturbation Theory rely on an accurate description of the Green functions 
beyond the linear regime. This is the purpose for instance of RPT and RegPT methods 
([28-30]). 
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